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Abstract 

We introduce quadri-tilings and show that they are in bijection with dimer models on a, family 
of graphs {R*} arising from rhombus tilings. Using two height functions, we interpret a sub- 
family of all quadri-tilings, called triangular quadri-tilings, as an interface model in dimension 2+ 
2. Assigning "critical" weights to edges of R* , we prove an explicit expression, only depending on 
the local geometry of the graph R* , for the minimal free energy per fundamental domain Gibbs 
measure; this solves a conjecture of We also show that when edges of R* are asymptotically 
far apart, the probability of their occurrence only depends on this set of edges. Finally, we give 
an expression for a Gibbs measure on the set of all triangular quadri-tilings whose marginals are 
the above Gibbs measures, and conjecture it to be that of minimal free energy per fundamental 
domain. 

1 Introduction 

In this paper, we introduce quadri-tilings and a sub-family of all quadri-tilings called triangular 
quadri-tilings. In order to explain the originality of this model, let us go back to the yet classical 
domino and lozenge tilings, see for example [S1E1IZ1E]- Both models are dimer models (see below) 
on a fixed graph, the square lattice Z 2 , and the equilateral triangular lattice T, respectively. By 
the means of a height function, they can be interpreted as random discrete interfaces in dimension 

2 + 1, that is as random discrete surfaces of dimension 2 in a space of dimension 3 that have been 
projected to the plane ^Sj. Keeping this in mind, one can now explain the interesting feature of 
quadri-tilings: they correspond to dimer models on a family of graphs, instead of a fixed graph as 
was the case up to now. Moreover, by the means of two height functions, triangular quadri-tilings 
can be interpreted as random interfaces in dimension 2 + 2, that is random surfaces of dimension 2 
in a space of dimension 4, that have been projected to the plane. It is the first such model arising 
from dimer models. Using tools of the dimer model, we are able to give an explicit expression for a 
Gibbs measure on triangular quadri-tilings, as well as a surprising property of the asymptotics of 
this measure. In the course of doing so, we prove a conjecture of |Hj. 

Quadri-tilings are defined as follows. Consider the set of right triangles whose hypotenuses have 
length two. Color the vertex at the right angle black, and the other two vertices white. A quadri- 
tile is a quadrilateral obtained from two such triangles in two different ways: either glue them 
along the hypotenuse, or supposing they have a leg of the same length, glue them along this edge 
matching the black (white) vertex to the black (white) one, see Figure ^ Note that both types of 
quadri-tiles have four vertices. A quadri-tiling of the plane is an edge-to-edge tiling of the plane 
by quadri-tiles that respects the coloring of the vertices, that is black (resp. white) vertices are 

*beatrice . detiliereSmath .unizh. ch 



1 



matched to black (resp. white) ones. An example of quadri-tiling is given in Figure ^ I n an that 
follows, we consider quadri-tilings of the plane that use finitely many different quadri-tiles up to 
isometry. 



The goal of Section |2] is to precisely describe the features of quadri-tilings. In order to give some 
insight, let us define 2-tiling models or equivalently dimer models. A 2-tile of an infinite graph 
G is a polygon made of two adjacent inner faces of G, and a 2-tiling of G is a covering of G with 
2-tiles, such that there are no holes and no overlaps. The dual graph G* of G is the graph whose 
vertices correspond to faces of G, two vertices of G* being joined by an edge if the corresponding 
faces are adjacent. A dimer configuration of G* , also called perfect matching, is a subset of 
edges of G* which covers each vertex exactly once. Then 2-tilings of the graph G are in bijection 
with dimer configurations of the dual graph G* , as explained by the following correspondence. 
Denote by /* the dual vertex of a face G, and consider an edge f*g* of G* . We say that the 2-tile 
of G made of the adjacent faces / and g is the 2-tile corresponding to the edge f*g*. Then, 
2-tiles corresponding to edges of a dimer configuration of G* form a 2-tiling of G. Let us denote 
by A4(G*) the set of perfect matching of the graph G* . 

Prior to describing Section[2J we need one more definition. If R is a rhombus tiling of the plane, then 
the corresponding rhombus- with-diagonals tiling, denoted by R, is the graph obtained from R 
by adding the diagonals of the rhombi. In the whole of this paper, we suppose that rhombus tilings 
of the plane have only finitely many rhombus angles. 

In Section 12 . 1| we prove the main feature of quadri-tilings, i.e. that they correspond to 2-tilings 
on a family of graphs, which consists of rhombus-with-diagonals tilings. More precisely, if T is a 
quadri-tiling, then by a geometric construction, we associate to T a rhombus-with-diagonals tiling 
R(T), such that T is a 2-tiling of R(T). The corresponding rhombus tiling R(T) is called the 
underlying rhombus tiling of T. 

Triangular quadri-tilings consist in the sub-family of all quadri-tilings whose underlying tiling is 
a lozenge tiling, where lozenges are defined to be 60°-rhombi; refer to Figure|2]for an example, and 
to Section 12.11 for the construction of the underlying tiling. In order to distinguish general rhombus 
tilings, denoted R, from lozenge tilings, we denote the latter by L. The set of all triangular quadri- 
tilings up to isometry is denoted by Q. Note that Q corresponds to two superposed dimer models. 




Figure 1: Two types of quadri-tiles (left), and a quadri-tiling (right). 
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Indeed, let T be a triangular quadri-tiling, then T is a 2-tiling of its underlying lozenge-with- 
diagonals tiling L(T), moreover the corresponding lozenge tiling L(T) is a 2-tiling of the equilateral 
triangular lattice T. 




Figure 2: Triangular quadri-tiling T (left), underlying lozenge- with-diagonals tiling L(T) (middle), 
corresponding lozenge tiling L(T) (right). 

Section 12.21 consists in the geometric interpretation of triangular quadri-tilings using height func- 
tions. On the vertices of every triangular quadri-tiling T, we define a Z- valued function h\, called 
the first height function, corresponding to the "height" of T interpreted as a 2-tiling of its 
underlying lozenge- with-diagonals tiling L(T). Then, we assign a second height function h,2 
(Thurston's height function on lozenges [T^]) corresponding to the height of L(T) interpreted as a 
2-tiling of T, see Figure HI Hence triangular quadri-tilings are characterized by two height func- 
tions, and so can be interpreted as discrete interfaces in dimension 2 + 2. In Section [2.31 we give 
elementary operations that allow to transform any triangular quadri-tiling of a simply connected 
region into any other. 

The dimer model belongs to the field of statistical mechanics, hence there are natural measures to 
consider, Boltzmann and Gibbs measures, which are defined as follows. Let G be an infinite graph, 
and let v be a positive weight function on the edges of G*. Consider a finite sub-graph G l of G, 
then the Boltzmann measure on the set of dimer configurations A4(G 1 *) of G 1 *, corresponding 
to the weight function u, is defined by 

am) = 

where Z(G 1 * , v) = Y^MeMlG 1 *) IleeM v i e ) ls * ne dimer partition function. A Gibbs measure 
is a probability measure on Ai(G*) with the following property. If the matching in an annular 
region is fixed, the matchings inside and outside of the annulus are independent of each other, and 
the probability of any interior matching M is proportional to f^eeA/ u ( e )- 

An important question in solving a dimer model is the study of local statistics, i.e. to obtain an 
explicit expression for the set of Gibbs measures. Kenyon, Okounkov and Sheffield give such 
an expression for the two-parameter family of Gibbs measures on dimer configurations of doubly 
periodic bipartite graphs. The expression they obtain involves the limiting inverse Kasteleyn ma- 
trix which is hard to evaluate in general, often implying elliptic integrals. In another paper [H], for 
graphs G which have bipartite duals and satisfy a geometric condition called isoradiality, Kenyon 
defines a specific weight function on the edges of G* called the critical weight function. He also 
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defines the Dirac operator K indexed by the vertices of G* , and gives an explicit expression for its 
inverse K~ l (see also Sections 13.11 and 13, 2(1 . The expression for K~ l has the interesting property 
of only depending on the local geometry of the graph. Kenyon conjectures that K" 1 is, in some 
sense to be determined, the limiting inverse Kasteleyn matrix. In Section |3J we consider a general 
rhombus-with-diagonals tiling of the plane R. It has the property of being an isoradial graph, so 
that we assign the critical weight function to edges of R* . Theorem 13.21 of Section 13.31 (see also 
Theorem 1 1 . 1 1 b elow) proves an explicit expression for a Gibbs measure fi R on A4(R*), as a function 
of K and K~ l . 

For every subset of edges e\ = W\b\, . . . , e^ = w^bk of R*, the cylinder {ex, . . . , e^} is defined to 
be the set of dimer configurations of R* that contain these edges. Then 

Theorem 1.1 There is a probability measure /j, r on A4(R*) such that, for every cylinders {ex, . . . , e&} 
ofR*, 

fj, R ( ei ,...,e k ) = (f[K(wi,bi) ] det {K-\b uWj )) . (1) 
Vf=i )i<*,3<k 

Moreover (j, R is a Gibbs measure on M(R*). When R* is doubly periodic, [i R is the unique Gibbs 
measure which has minimal free energy per fundamental domain among the two-parameter family 
of ergodic Gibbs measures of Jjj^. 

• Note that we do not ask the graph R* to be periodic. The proof of Theorem ll.ll is the subject 
of Section 0] The argument in the case where R* is not periodic relies on the argument in the 
case where R* is doubly periodic, combined with a non-trivial geometric property of rhombus 
tilings proved in Proposition ^. II "every finite simply connected sub-graph of a rhombus tiling 
can be embedded in a periodic rhombus tiling of the plane." 

• From the proof of Theorem II .!( it appears that the statement is true for all doubly periodic 
isoradial graphs with bipartite duals. Hence, Theorem 11.11 solves the conjecture of jS] of 
interpreting the inverse Dirac operator as the limiting inverse Kasteleyn matrix. Moreover, 
the fact that the measure fi R is of minimal free energy per fundamental domain makes it of 
special interest among the two-parameter family of ergodic Gibbs measures of [2]. 

• Using the locality property of K^ 1 mentioned above, we deduce that the expression Q only 
depends on the local geometry of the graph, hence it yields an easy way of computing local 
statistics explicitly. This is very surprising in regards of the expression obtained in [5], and 
we believe this locality property to be true only in the isoradial case with critical weights. 

In Section we extend the notion of Gibbs measure to the set of all triangular quadri-tilings Q. 
Then, as a corollary to Theorem II .11 we deduce an explicit expression for such a Gibbs measure \i, 
and conjecture it to be that of minimal free energy per fundamental domain among a four param- 
eter family of Gibbs measures. 

In Section we consider a general rhombus-with-diagonals tiling of the plane R. We assign the 
critical weight function to edges of R*, and let K be the Dirac operator indexed by vertices of 
R*. Theorem 16.11 of Section 16.11 fsee also Theorem 11.21 below) establishes that asymptotically (as 
\b — w\ — ► oo) and up to the second order term, K _1 {b, w) only depends on the rhombi to which 
the vertices b and w belong, and else is independent of the structure of the graph R. For a general 
isoradial graph, Kenyon jS] gives an asymptotic formula for K (b, w) which depends on the angles 
of an edge-path from w to b. Hence, it is an interesting and surprising fact that the dependence on 
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the edges along the path should asymptotically disappear in the case of rhombus-with-diagonals 
tilings. 

Theorem 1.2 As \b — w\ ^ oo, iT -1 (&, w) is equal to 

J_ ( 1 , e ~'(«l+«2) \ J_ / e 2ia 1+e 2ifl 2 e -.(39 1 +D 2 ) +e - i(ei +38 2 ) \ / 1 \ 

2tt \b-w "T 6-tjj ,1 + 2^ (6-iu) 3 r (6-iB) 3 / V I 6 "™! 3 / ' 

As a consequence of Theorem ll.21 we deduce that when edges ex, . . . , of R* are asymptotically far 
apart, /i R (ei, . . . , e&) only depends on the rhombi to which the edges e\, . . . , belong, and else is 
independent of the structure of the graph R (Corollary I6.2|h We conclude by giving a consequence 
of Corollary 16.21 for the measure (i on triangular quadri-tilings ( Corollary 16. 4j) . 

Acknowledgments: We thank Richard Kenyon for proposing the quadri-tiling model and asking the 
questions related to it; we are grateful to him for the many enlightening discussions. We also thank 
the referee for the many pertinent remarks which have helped to increase the quality of this paper. 

2 Features of quadri-tilings 

2.1 Underlying rhombus-with-diagonals tilings 

Lemma 2.1 Quadri-tilings are in one-to-one correspondence with 2-tilings of graphs which are 
rhombus-with-diagonals tilings of the plane. 

Proof: Consider a quadri-tiling of the plane T. Denote by R the tiling of the plane obtained 
from T by drawing, for each quadri-tile, the edge separating the two right triangles. Let b be a 
black vertex of R, denote by w\, . . . , the neighbors of b in cclw (counterclockwise) order. In each 
right triangle, the black vertex is adjacent to two white vertices, and since the gluing respects the 
coloring of the vertices, w±, . . . ,Wk are white vertices. Moreover, b is at the right angle, so k = 4 
and the edges W1W2, W2W3, w^w± are hypotenuses of right triangles. Therefore W\, . . . ,104 

form a side-length-2 rhombus, and b stands at the crossing of its diagonals. This is true for any 
black vertex b of R, so R is a rhombus-with-diagonals tiling of the plane, and T is a 2-tiling of R. 

□ 

As a consequence of Lemma 12.11 a quadri-tiling T is a 2-tiling of a unique rhombus-with-diagonals 
tiling, which we call the underlying rhombus-with-diagonals tiling, and denote by R(T), see 
Figure [21 

2.2 Height functions 

We define a first height function h\ on vertices of every quadri-tiling T. Moreover, when T is 
a triangular quadri-tiling, we define a second height function h,2 on vertices of T. Using h\ 
and h,2, we interpret triangular quadri-tilings as discrete 2-dimensional surfaces in a 4-dimensional 
space projected to the plane. 

2.2.1 First height function 

Consider a quadri-tiling of the plane T, then T is a 2-tiling of its underlying rhombus-with-diagonals 
tiling R(T). In order to define the first height function hi, we need a bipartite coloring of the faces 
of R(T), which is given by the following. 
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Lemma 2.2 Let R be a rhombus tiling of the plane, and R be the corresponding rhombus-with- 
diagonals tiling. Then R has a bipartite coloring of its faces which is also a bipartite coloring of 
the vertices of R* . 

Proof: Cycles corresponding to the faces of the graph R have length four, thus R has a bipartite 
coloring of its vertices, say black and white. Consider a face of R and orient its boundary edges 
cclw. If the white vertex of the hypotenuse-edge comes before the black one, assign color black to 
the face, else assign color white. This defines a bipartite coloring of the faces of R, which is also a 
bipartite coloring of the vertices of R* (see Figure |2J) . 



Figure 3: Bipartite coloring of the vertices of R (left), and corresponding bipartite coloring of the 
faces of R and of the vertices of R* (right). 



Consider the bipartite coloring of the faces of R(T). Orient the edges around the black faces cclw, 
edges around the white faces are then oriented cw, and define hi on the vertices of T as follows. Fix 
a vertex v\ on a boundary edge of a rhombus of R(T), and set h\(v\) = 0. For every other vertex v 
of T, take an edge-path 71 from v\ to v which follows the boundaries of the quadri-tiles of T. The 
first height function hi changes by ±1 along each edge of 71: if an edge is oriented in the direction 
of the path, then hi increases by 1, if it is oriented in the opposite direction, then hi decreases by 
1. The value hi(v) is independent of the path 71 because the plane is simply connected, and the 
height change around any quadri-tile is zero. An example of computation of hi is given in Figure 



The following lemma gives a bijection between 2-tilings of a rhombus-with-diagonals tiling R and 
first height functions defined on vertices of R. 

Lemma 2.3 Fix a vertex v\ on a boundary edge of a rhombus of R. Let hi be a "L-valued function 
on the vertices of R satisfying the following two conditions: 

• h\{vi) = 0, 

• hi(v) = hi(u) + 1, or hi(v) = hi(u) — 2, for any edge uv oriented from u to v. 

Then, there is a bijection between functions hi satisfying these two conditions and 2-tilings of R. 

Proof: The idea of the proof closely follows 3 j. If T is a 2-tiling of R, then the first height 
function defined above satisfies the two conditions of the lemma: if an edge uv, oriented from u 
to v, belongs to the boundary of a quadri-tile, it satisfies hi{v) = hi(u) + 1, else if it lies across a 
quadri-tile, it satisfies hi(v) = hi(u) — 2. 

Conversely, consider a Z-valued function hi as in the lemma. Then, anytime there is an edge uv 



□ 
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satisfying \hi(v) — h\{u)\ = 2, put a quadri-tile made of the two right triangles adjacent to this 
edge. This defines a 2-tiling of R. □ 



2.2.2 Second height function 

Consider a triangular quadri-tiling T. Let L{T) be its underlying lozenge- with-diagonals tiling, and 
h\ be the first height function on vertices of T. The lozenge tiling L(T) corresponding to L(T) is 
a 2-tiling of the equilateral triangular lattice T. Moreover T has a bipartite coloring of its faces, 
say black and white. Orient the edges around the black faces cclw, edges around the white faces 
are then oriented cw. Thurston |18j defines the second height function hi as follows: choose a 
vertex vi of L(T), and set hi(vi) = 0. For every other vertex v of L(T), take an edge-path 72 
from V2 to v which follows the boundaries of the lozenges of L(T). The second height function hi 
changes by ±1 along each edge of 72: if an edge is oriented in the direction of the path, then hi 
increases by 1, if it is oriented in the opposite direction, then hi decreases by 1. The value hi(v) is 
independent of the path 72. For convenience, we choose vi to be the same vertex as v±, and denote 
this common vertex by t>o, so that hi(vo) = hi(vo) = 0. An analog to Lemma EH1 gives a bijection 
between second height functions and lozenge tilings of the plane, hence we deduce that triangular 
quadri-tilings are characterized by hi and hi. 

Let us define a natural value for the second height function at the vertex in the center of the 
lozenges of L(T). When going cclw around the vertices of a lozenge t of L(T), starting from the 
smallest value of hi say h, vertices take on successive values h, h + 1, h + 2, h + 1, so that we assign 
value h + 1 to the vertex in the center of the lozenge i. An example of computation of hi is given 
in Figure 



1 
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Figure 4: Triangular quadri-tiling with height functions h\ (above) and hi (below). 

In Thurston's geometric interpretation ^H], a lozenge tiling is seen as a surface S in Z 3 (where 
the diagonals of the cubes are orthogonal to the plane) that has been projected orthogonally to 
the plane. The surface S is determined by the height function hi. In a similar way, a triangular 
quadri-tiling of the plane T can be seen as a surface S\ in a 4-dimensional space that has been 
projected to the plane; S\ can also be projected to Z 3 (Z 3 is the space Z 3 where cubes are drawn 
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with diagonals on their faces), and one obtains a surface S2. When projected to the plane, S2 is 
the underlying lozenge-with-diagonals tiling L(T). 

2.3 Elementary operations 

Consider a finite simply connected sub- graph G of the equilateral triangular lattice T, and let dG 
be the cycle of G consisting of its boundary edges. Denote by Q(dG) the set of triangular quadri- 
tilings whose underlying tilings are lozenge tilings of G. Let L 1 be a lozenge tiling of G, and L 1 be 
the corresponding lozenge-with-diagonals tiling. Then using the bijection between the first height 
function and 2-tilings of L 1 we obtain, in exactly the same way as Elkies, Kuperberg, Larsen, Propp 
[3] have for domino tilings, the following lemma: 



Lemma 2.4 Every 2-tiling of L can be transformed into any other by a finite sequence of the 
following operations, (in brackets is the number of possible orientations for the graph corresponding 
to the operation): 




Figure 5: Quadri-tile operations. 



Let us call quadri-tile operations the 21 operations described in Lemma 12.41 Moreover, every 
lozenge tiling of G can be transformed into any other by a finite sequence of lozenge operations: 





Figure 6: Lozenge operations. 
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Note that if L is any lozenge tiling of G, then L is quadri-tilable with quadri-tiles obtained by 
cutting in two every lozenge along one of its diagonals. Moreover, when one performs a lozenge 
operation on such a quadri-tiling, one still obtains a quadri-tiling of Q{dG). Let us call elementary 
operations the quadri-tile operations and the lozenge operations performed on quadri-tilings as 
described above. Then we have: 

Lemma 2.5 Every quadri-tiling of Q(dG) can be transformed into any other by a finite sequence 
of elementary operations. 

Proof: This results from Lemma l2.4| and the above observation. □ 



3 Gibbs measure on quadri-tilings 

This section aims at giving a precise statement of Theorem 11.11 of the introduction (see Theorem 
13.21 below). Sections 13.11 and 13.21 are taken form [S] and give a precise definition of an isoradial 
graph, the critical weight function, the Dirac and inverse Dirac operator. Section 13.31 consists in 
the statement of Theorem 13.21 

3.1 Isoradial graphs and critical weight function 

The definition of the critical weight function follows jS] . It is defined on edges of graphs satisfying 
a geometric condition called isoradiality: all faces of an isoradial graph are inscribable in a circle, 
and all circumcircles have the same radius. 

Note that if R is a rhombus tiling of the plane, then the corresponding rhombus-with-diagonals 
tiling R is an isoradial graph. Let us consider the embedding of the dual graph R* (the same 
notation is used for the one-skeleton of a graph and its embedding) where the dual vertices are the 
circumcenters of the corresponding faces. Then R* is also an isoradial graph and the circumcenters 
of the faces are the vertices of R. 

To each edge e of R*, we associate a unit side- length rhombus R(e) whose vertices are the vertices 
of e and the vertices of its dual edge. Let R = U eg i?*i?(e). Note that the dual edges corresponding 
to the boundary edges of the rhombi of R have length zero, and that the rhombi associated to these 
edges are degenerated. 

For each edge e of R*, define u{e) = 2sin#, where 29 is the angle of the rhombus R(e) at the vertex 
it has in common with e; 9 is called the rhombus angle of the edge e. Note that u(e) is the length 
of e*, the dual edge of e. The function v is called the critical weight function. 

3.2 Dirac and inverse Dirac operator 

Results and definitions of this section are due to Kenyon [B] , see also Mercat Jl] . Note that this 
section (as the previous one) is true for general isoradial graphs with bipartite dual graphs. 

Let R be a rhombus tiling of the plane, then by Lemma 12.21 R* is a bipartite graph. Let B (resp. 
W) be the set of black (resp. white) vertices of R*, and denote by v the critical weight function 
on the edges of R* . The Hermitian matrix K indexed by the vertices of R* is defined as follows. 
If vi and V2 are not adjacent K(p\,v%) = 0. If w £ W and b G B are adjacent vertices, then 
K(w, b) = K(b, w) is the complex number of modulus v{wb) and direction pointing from w to b. 
If w and b have the same image in the plane, then \K(w,b)\ = 2, and the direction of K(w,b) is 
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that which is perpendicular to the corresponding dual edge, and has sign determined by the local 
orientation. The infinite matrix K defines the Dirac operator K : C V ( R *' — > C V ( R *\ by 

(Kf){v)= Y, K(v,u)f(u). 
ueB* 

where V{R*) denotes the set of vertices of the graph R*. 

The inverse Dirac operator K~ l is defined to be the operator satisfying: 

1. KK- 1 = Id, 

2. K~ l (b, w) — > 0, when \b — w\ — > oo. 

Kenyon [S] obtains an explicit expression for K^ 1 . Before stating his theorem, we need to define 
the rational functions f wv (z). Let w be a white vertex of R*. For every other vertex v, define 
f wv (z) as follows. Let w = vq, vi, i?2, • • • , Ufe = v be an edge-path of R, from w to v. Each edge 
VjVj+x has exactly one vertex of R* (the other is a vertex of R). Direct the edge away from this 
vertex if it is white, and towards this vertex if it is black. Let e taj be the corresponding vector in 
R (which may point contrary to the direction of the path), then f wv is defined inductively along 
the path, starting from 

fww{.z) — L 

If the edge leads away from a white vertex, or towards a black vertex, then 

j, , n _ fwvj (z) 

Jwv ]+1 {Z) — - _ ^ ia , , 

else, if it leads towards a white vertex, or away from a black vertex, then 

fwvj + ii z ) = fwvj{z).(z — e 3 '). 

The function f wv (z) is well defined (i.e. independent of the edge-path of R from w to v), because 
the multipliers for a path around a rhombus of R come out to 1. For a black vertex b the value 
K~ l {b, w) will be the sum over the poles of f w b(z) of the residue of /„& times the angle of z at the 
pole. However, there is an ambiguity in the choice of angle, which is only defined up to a multiple 
of 2-7T. To make this definition precise, angles are assigned to the poles of f w b(z). Working on the 
branched cover of the plane, branched over w, so that for each black vertex b in this cover, a real 
angle #o is assigned to the complex vector b — w, which increases by 2tt when b winds once around 
w. In the branched cover of the plane, a real angle in [8q — tt + A, #o + it — A] can be assigned to 
each pole of f w b, for some small A > 0. 

Theorem 3.1 jH] There exists a unique K~ l satisfying the above two properties, and K" 1 is given 
by: 

K~ l (b, w) = J f w b(z) log z dz, 

where C is a closed contour surrounding cclw the part of the circle {e ld \9 € [9q — ir + A, 9$ + tt — A]}, 
which contains all the poles of f w b, and with the origin in its exterior, see Figure^ 

The remarkable property of K~ l {b, w) is that it only depends on the local geometry of the graph, 
i.e. on an edge-path from w to 6. 
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e 




Figure 7: An example of contour C. 

3.3 Statement of result 

Let R be a rhombus tiling of the plane, and R be the corresponding rhombus-with-diagonals tiling. 
Suppose that the critical weight function v is assigned to edges of R* , and let K be the Dirac 
operator indexed by the vertices of R*. Moreover, recall that if e% = w\b\, . . . , = w^bk is a subset 
of edges of R* , then the cylinder {ei, . . . , e^} is defined to be the set of dimer configurations of R* 
which contain these edges. Let A be the field consisting of the empty set and of the finite disjoint 
unions of cylinders. Denote by a(A) the cr-field generated by A. 

Theorem 3.2 There is a probability measure on (A4(R*),a(A)) such that for every cylinder 
{ei,...,e fc } of R*, 

fi R (ei,...,e k ) = njK(wi,bi)J ^det^K^ib^Wj)) . (2) 

Moreover [i R is a Gibbs measure on A4(R*). When R* is doubly periodic, fi R is the unique Gibbs 
measure which has minimal free energy per fundamental domain among the two-parameter family 
of ergodic Gibbs measures of J3j. 

From the proof, it appears that Theorem 13.21 is in fact true for all doubly periodic isoradial graphs 
with bipartite dual graphs. We refer to the introduction for comments on Theorem 13.21 



4 Proof of Theorem EL2] 

The proof of Theorem 13. 21 uses Propositions 14. ll and !4.9l below. Proposition 14. II is a geometric prop- 
erty of rhombus tilings and is the subject of Section f4. II Proposition 14.91 concerns the convergence 
of the Boltzmann measure on some appropriate toroidal graphs, it is the subject of Section 14.31 
In Section 14.21 we introduce the real Dirac operator and its inverse. This operator is related to 
the Dirac operator and is needed for the proof of Proposition 14.91 Theorem 13.21 is then proved in 
Section |OJ 
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4.1 Geometric property of rhombus tilings 

Proposition 4.1 Let R be a rhombus tiling of the plane, then any finite simply connected sub-graph 
P o/R can be embedded in a periodic rhombus tiling S of the plane. 

Proof: This proposition is a direct consequence of Lemmas 14.21 14.31 and Theorem 14 . 41 b elow . □ 

The notion of train-track has been introduced by Mercat JH]) see also Kenyon and Schlenker 
[51 111 j . A train-track of a rhombus tiling is a path of rhombi (each rhombus being adjacent along 
an edge to the previous rhombus) which does not turn: on entering a rhombus, it exits across 
the opposite edge. Train-tracks are assumed to be maximal in the sense that they extend in both 
directions as far as possible. Thus train-tracks of rhombus tilings of the plane are bi-infinite. Each 
rhombus in a train-track has an edge parallel to a fixed unit vector e, called the transversal 
direction of the train-track. Let us denote by t e the train-track of transversal direction e. In 
an oriented train-track (i.e. the edges of the two parallel boundary paths of the train-track have 
the same given orientation), we choose the direction of e so that when the train-track runs in the 
direction given by the orientation, e points from the right to the left. The vector e is called the 
oriented transversal direction of the oriented train-track. A train-track cannot cross itself, 
and two different train-tracks can cross at most once. A finite simply connected sub-graph P of a 
rhombus tiling of the plane R is train-track-convex, if every train-track of R that intersects P 
crosses the boundary of P twice exactly. 

Lemma 4.2 Let R be a rhombus tiling of the plane, then any finite simply connected sub-graph P 
of R can be completed by a finite number of rhombi of R in order to become train-track- convex. 

Proof: Let e±, . . . , e m be the boundary edges of P. Every rhombus of P belongs to two train- 
tracks of R, each of which can be continued in both directions up to the boundary of P. In both 
directions the intersection of each of the train-tracks and the boundary of P is an edge parallel to 
the transversal direction of the train-track. Thus, to take into account all train-tracks of R that 
intersect P, it suffices to consider for every i the train-tracks t ei associated to the boundary edges 
of P. Consider the following algorithm (see Figure [SJ. 

Set Qi = P. 

For i = 1, . . . , m, do the following: 

Consider the train-track t ei , and let 2n« be the number of times t ei intersects the bound- 
ary of Qi. 

• If n>i > 1: there are rtj — 1 portions of t ei that are outside of Qj, denote them by 
t ei , . . . ,t ei ni ~ l . Then, since Qi is simply connected, for every j, R\ (Qj U t e J) is 
made of two disjoint sub-graphs of R, one of which is finite (it might be empty 
in the case where one of the two parallel boundary paths of t ei 3 is part of the 
boundary path of Qj). Denote by g ei 3 the simply connected sub-graph of R made 
of the finite sub-graph of R \ (Qi U t e J) and of t ei 3 . Denote by b e 3 the portion of 
the boundary of Qj which bounds g ei 3 . Replace Qj by Qj+i = Qj U (U™!^ 1 ^/). 
By this construction t ei intersects the boundary of Qj+i exactly twice, and Qj+i is 
simply connected. 
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• If m = 1: set Qj+i = Qi. 




Figure 8: One step of the algorithm. 

Let us show that at every step the train-tracks of R that intersect Qj and Qj+i are the same. By 
construction, boundary edges of Qj+i are boundary edges of Qj and of t e J , for every j. Let / be 
an edge on the boundary of Qj+i, but not of Qj, that is / is on the boundary of t ei 3 for some j, 
thus tf crosses g e J . Since two train-tracks cross at most once, tf has to intersect b ei 3 , which means 
tf also crosses Qj. From this we also conclude that if a train-track intersects the boundary of Qj 
twice, then it also intersects the boundary of Qj + i twice. 

Thus all train-tracks that intersect Q m +i cross its boundary exactly twice, and Q m +i contains P. □ 



Lemma 4.3 Let R be a rhombus tiling of the plane. Then any finite simply connected train-track- 
convex sub-graph P o/R can be completed by a finite number of rhombi in order to become a convex 
polygon Q, whose opposite boundary edges are parallel. 

Proof: Let e\, . . . , e m be the boundary edges of P oriented cclw. Since P is train-track-convex, 
the train-tracks t ei , . . . , t em intersect the boundary of P twice, so that there are pairs of parallel 
boundary edges. Let us assume that the transversal directions of the train-tracks are all distinct 
(if this is not the case, one can always perturb the graph a little so that it happens). Let us also 
denote by t ei , . . . , t em the portions of the bi- infinite train-tracks of R in P. In what follows, indices 
will be denoted cyclically, that is ej = e m+ j. Write Xj (resp. i/j) for the initial (resp. end) vertex 
of an edge ej. 

Let ej,ej + i be two adjacent boundary edges of P. Consider the translate e\ +1 of ej+i so that the 
initial vertex of e* +1 is adjacent to the initial vertex of ej. Then we define the turning angle from 

ej to ej+i (also called exterior angle) to be the angle eje* +1 , and we denote it by eijCi+1 . If ej,ej 
are two boundary edges, then the turning angle from ej to ej is defined by Y?a=i @e a ,e a+1 , and 
is denoted by <9 ei>ei . 

Properties 

L Ea=l 8e a ,e a+1 = 2vr. 
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2. If ej, ej are two boundary edges, and if 7 = . . . , / n } is an oriented edge-path in P from 
Vi to Xj , then 6 e ^ ej = euh + ^=1 + %n,e r 

3. If e« is a boundary edge of P, and is the second boundary edge at which t ei intersects the 
boundary of P, then ei ,e k = 7T- Thus and ej are oriented in the opposite direction, and we 
denote by ef 1 . 

4. P is convex, if and only if every train-track of P crosses every other train-track of P. 

We first end the proof of Lemma l4.3| and then prove Properties 1 to 4. 

Note that Properties 1 and 2 are true for any finite simply connected sub- graph of R. 

The number of train-tracks intersecting P is n = m/2. So that if every train-track crosses every 

other train-track, the total number of crossings is n(n — l)/2. Consider the following algorithm 

(see Figure El for an example). 

Set Qi = P, n\ = the number of train-tracks that cross in Qi. 
For i = 1,2,... do the following: 

• If rii = n(n — l)/2: then by Property 4, Qj is convex. 

• If i%i < n(n — l)/2: then by Property 4, e -.,e-.+i < for some j% £ {!>••■ ; m }- 
Add the rhombus £j t of parallel directions e^, ej i+ \ along the boundary of Qj. Set 
Qi_l_i = Qj U £j i , and rename the boundary edges e\, . . . , e m in cclw order. Then 
the number of train-tracks that cross in Qj+i is n% + 1, set rij + i = nj + 1. Note 
that if Property 4 is true for Qj, it stays true for Qj+i, and note that the same 
train-tracks intersect Qj and Qj+i. 




Figure 9: Example of application of the algorithm. 

For the algorithm to be able to add the rhombus £j. at every step, we need to check that: 

V i, 9 e . , e . , , > — 7T, and 9 e . e . , „ > — tt. (3) 

Assume we have proved that for any finite simply connected train-track-convex sub-graph P of R 
we have: 

Vt,i, f ei , e , > -7T. (4) 
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Then Properties 1 and 2 imply that if (jlj) is true for Qj it stays true for Qj+i, moreover @ implies 
(jSJ). So let us prove (J3J by induction on the number of rhombi contained in P. If P is a rhombus, 
then Q is clear. Now assume P is made of k rhombi. Consider the train-tracks in P adjacent to 
the boundary (every boundary edge e of P belongs to a rhombus of P which has parallel directions 
e and /; for every boundary edge e, the train-track of transversal direction / is the train-track 
adjacent to the boundary). Denote the train-tracks adjacent to the boundary by t\, . . . ,t p in cclw 
order, and write fp for the oriented transversal direction of tp (when the boundary edge-path of 
P is oriented cclw). Consider two adjacent boundary edges ei,ei + i of P that don't belong to the 
same boundary train-track. That is e« belongs to tp, and e^+i to t/3+1. Then either fpfp+i < 

or fpfp+i > 0, in the second case tp and tp+i cross and their intersection is a rhombus ip of P. 
The rhombus ip has boundary edges ej,ej + i, and fp+i = e{~ 1 ,fp~ 1 = ej + i -1 . Now Property 1 

implies that X^/3=i fpfp+i = 2tt, so that there always exists (3q such that fp fp +\ > 0. Removing 
£p from P and using the assumption that P is train-track-convex, we obtain a graph P' made of 
k — 1 rhombi which is train-track-convex. By induction, 8 e j > —it for every boundary edge of P', 
and using Property 2, we conclude that this stays true for P. 

Denote by Q the convex polygon obtained from P by the algorithm, and assume that opposite 
boundary edges are not parallel. Then there are indices i and j such that comes before ej, 
and ej 1 comes before e," 1 . This implies that 6 eu£j = —0 e .-i e .-i, so that one of the two angles is 
negative, which means Q can not be convex. Thus we have a contradiction, and we conclude that 
opposite boundary edges of Q are parallel. 

Proof of Properties 1 to 4- 
1. and 2. are straightforward. 

3. When computing # ei ,e fe along the boundary edge-path of t £i we obtain ir, so by Property 2 we 

deduce that ei ,e k = vr in P. 

4. P is convex if and only if, for every i, 9 ei ,e i+ i > 0, which is equivalent to saying that, for every 

i ^ j, e ^ e . > 0. Therefore Property 4 is equivalent to proving that 9 euej > 0, for every i ^ j, 
if and only if every train-track of P crosses every other train-track of P. 
Assume there are two distinct train-tracks t ee and t ek that don't cross in P. Then, in cclw 
order around the boundary of P, we have either ee, e^ -1 , e^, ef \ or eg, e^, e^ -1 , ef x . It 
suffices to solve the second case, the first case being similar. By Property 1, Q ei ,e k +$e fe ,e fc - 1 + 
8e k -\e r ^ + Qe t -\e t = 271 *■ Moreover by Property 3, 6> efc )6fc -i = et -i )Bjt = it, which implies 
&e e ,e k = ~&e k ~ 1 ,e e - 1 - Since all train-tracks have different transversal directions, either ei)( , k 
or G ek -i^ ee -i is negative. 

Now take two boundary edges ei,ej of P (with i ^ j, and ej ^ e^ -1 ), and assume the 
train-tracks t ei ,t e . cross inside P. Then in cclw order around the boundary of P, we have 
either ej, ej , e^ -1 , ej, or ej, ej, ej -1 , ey . It suffices to solve the second case since the first 
case can be deduced from the second one. The intersection of t e . and t e . is a rhombus t. Let 
ej 1 (resp. e~ ) be the boundary edge of £ parallel and closest to ej (resp. e^), oriented in 
the opposite direction, then 9~-i ~-i < 0. Let jj (resp. ji) be the boundary edge-path of 

t e . (resp. t ei ) from yj to Xj (resp. from jji to Xj), and let Q be the sub-graph of R whose 
boundary is e^, ej+i, . . . ,ej,jj, ej , e" 1 , 7^. Since t ei and t ej intersect the boundary of P twice, 
they also intersect the boundary of Q twice. Moreover t e . and t e . don't cross in Q, so that 
= -6L_! -1 > 0. □ 

□ 
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Theorem 4.4 Any convex 2n-gon Q whose opposite boundary edges are parallel and of the same 
length can be embedded in a periodic tiling of the plane by Q and rhombi. 

Proof: Let e±, . . . , e n , e\~ , . . . , e n _1 be the boundary edges of the polygon Q oriented cclw. If 
n < 3, then Q is either a rhombus or a hexagon, and it is straightforward that the plane can be 
tiled periodically with Q. 

If n > 4, for k = 1, . . . , n — 3, do the following (see Figure ITT)]): along e n -t add the finite train-track 
t en _ k of transversal direction e n _fc, going away from Q, whose boundary edges starting from the 
boundary of Q are: 

ei,e2,ei,e 3 ,e 2 ,ei,. . . , e ra _ fc _ 2 , • • • , e^ . 




Figure 10: Fundamental domain of a periodic tiling of the plane by dodecagons and rhombi. 

Since the polygon Q is convex, the rhombi that are added are well defined, moreover the intersection 
of t ei and the boundary of Q is the edge e^, and t ei doesn't cross t e . when i ^ j. So we obtain 
a new polygon Q' made of Q and rhombi, whose boundary edge-path is 71, ... ,76, (when starting 
from the edge e n _1 of Q), where: 

71 = en -1 , ex, e 2 ,jzx, e 3 , e 2 , ei , . . . , e n _ 3 , . . .,ey 

72 = e„- 2 , • • • , ei, 

7 3 = e n _i, 

74 = ei" 1 , . . . , e n _ 3 ~\ . . . , ei~\e 2 ~V i~ 1 ,e w , 

75 = ei- 1 ,...,e n _ 2 _1 , 

76 = e n _i _1 . 

Noting that 74 = 7i _1 ,75 = 7 2 ~ 1 >76 = 73 _1 ) and using the fact that the plane can be tiled period- 
ically with hexagons which have parallel opposite boundary edges, we deduce that the plane can 
be tiled with Q', that is it can be tiled periodically by Q and rhombi. □ 
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Remark 4.5 It was pointed out by the referee that Theorem \4-4\ might be known already. After 
a second look at the extensive literature on tilings, we were not able to find a reference, but any 
information is of course welcome. 

4.2 Real Dirac and inverse real Dirac operator 

In the whole of this section we let R be a rhombus tiling of the plane, and R be the corresponding 
rhombus-with-diagonals tiling. Assume that the critical weight function v is assigned to edges of 
R*, and denote by K the Dirac operator indexed by the vertices of R*. 

The proof of Theorem \'A.'2\ requires to take the limit of Boltzmann measures on some appropriate 
toroidal graphs (see Section r4.3|) . In order to do this, we need to introduce the real Dirac operator 
denoted by K. Both the Dirac operator K and the real Dirac operator K are represented by infinite 
weighted adjacency matrices indexed by the vertices of R*. For K, the edges of R* are un-oriented 
and weighted by their critical weight times a complex number of modulus 1. For K, edges of R* are 
oriented with a clockwise odd orientation, and are weighted by their critical weight. Both weight 
functions yield the same probability measure on finite simply connected sub-graphs of R*, but, and 
this is the reason why we introduce the real Dirac operator, these weights do not yield the same 
probability distribution on toroidal sub-graphs of R*. 

The structure of this section is close to that of Section fH. 21 We first define the real Dirac operator 
K and its inverse K _1 and then, using results of [S], we prove the existence and uniqueness of K _1 
by giving an explicit expression for K _1 . Note that this section is actually true for general isoradial 
graphs with bipartite dual graphs. 

4.2.1 Real Dirac operator 

Let us define an orientation of the edges of R*. An elementary cycle C of R* is said to be clockwise 
odd if, when traveling cw around the edges of C, the number of co-oriented edges is odd. Note 
that since R* is bipartite, the number of contra-oriented edges is also odd. Kasteleyn jS] defines 
the orientation of the graph R* to be clockwise odd if all elementary cycles are clockwise odd. 
He also proves that, for planar simply connected graphs, such an orientation always exists. 
Consider a clockwise odd orientation of the edges of R* . Define K to be the infinite adjacency matrix 
of the graph R*, weighted by the critical weight function v. That is, if v\ and V2 are not adjacent, 
K(wi,w 2 ) = 0. Hw G W and b G B are adjacent vertices, then K(w,b) = —K(b,w) = (—l) l ( w ' bS >v(wb), 
where h w ,b) = if the edge wb is oriented from w to 6, and 1 if it is oriented from b to w. The 
infinite matrix K defines the real Dirac operator K: C V ( R *^ — > C^^***, by 

("</)(«)= £ K(u, «)/(«), 

The matrix K is also called a Kasteleyn matrix for the underlying dimer model. 

4.2.2 Inverse real Dirac operator 

The inverse real Dirac operator K _1 is defined to be the unique operator satisfying 
1. KK- 1 = Id, 
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2. K _1 (6, w) -> 0, when \b 



W — > OO. 



Let us define the rational functions f wx (z). They are the analogous of the rational functions f wv (z), 
but are defined for vertices x S R* (whereas the functions f wv (z) were defined for vertices v € R). 
Let w € W, and let x € -B (resp. a; € W); consider the edge-path w = wi, b\, . . . , Wk, bk = x (resp. 
w = wi, bi, ... , tUfc, bk,Wk+i = x) of i?* from w to x. Let R(wjbj) be the rhombus associated to the 
edge Wjbj, and denote by Wj,Xj, bj,yj its vertices in cclw order; e m i is the complex vector yj — Wj, 
and e^' is the complex vector xj — Wj. In a similar way, denote by Wj+i,x'j,bj,y'j the vertices of 

the rhombus R(wj + \bj) in cclw order, then e a i is the complex vector y'- — Wj+i, and e*^ is the 
complex vector a;'- — tOj+i. Then f U)a: (^) is defined inductively along the path, 

^ww \Z) 1) 

W(*) = f ^( z )(— to^jr^y 

W +1 (z) = ^(^(-l^+i-^e-^Cz-e^JCz-e^). 



Remark 4.6 

VFe /iaue i/ie following relations between the real and the complex case. 
1. VweW,VxeBuW, f wx (z) = W0)/«(4 

£ Viu £ W, VfeG -B, smc/i i/iaf u; is adjacent to b, K(w,b) = f w b(0)K(w,b). 
Proof: 

1. This is a direct consequence of the definitions of the functions f wx and f wx . 

2. Let R(wb) be the rhombus associated to the edge wb, and let w, x, b, y be its vertices in cclw 
order. Denote by e ia the complex vector y — w, and by e 1 ^ the complex vector x — w. Let be the 
rhombus angle of the edge wb. By definition we have, 

i O.—0 - a— /3 

K(w,b) = (-1) W) 2 sin 9 = (-l)W) e 2 ~ e L 

tf(u>,6) = i(e 4/3 -e ia ), 
WO) = (-l)W) e -^. 

Combining the above three equations yields 2. □ 



Lemma 4.7 TTie function f wx is well defined. 

Proof: Showing that the function f wx is well defined amounts to proving that f wx is inde- 
pendent of the edge-path of R* from w to x. This is equivalent to proving the following: let 
wi, bi, . . . , Wk, bk, Wk+i = w\ be the vertices of an elementary cycle C of R* , where vertices are 
enumerated in cclw order; if f WlWl (z) = 1 then f Wl w k+1 (z) = 1. Let us use the notations introduced 
in the definition of f wx , and denote indices cyclically, that is k + 1 = 1. By Remark 14.61 we have 



'w 1 w k+1 \ z ) — >wiw k+1 (UJ Jwiw k+1 
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Since the function f wx is well defined, f wlWh+1 (z) = 1. Hence, it remains to prove that f W iw k+1 {ty 
1 . By definition of f wx , we have 

W +1 (0) = n ((-l) (I( ^'^ )+I( ^+ 1 '^ )) e l '^e-^] . 

Moreover for every j, a'j = 0j, so that 

f OT+1 (0) = n ((-l) {I ^^ )+ Wi^)) e ^ e ^') . 

Let 0j (resp. 0^) be the rhombus angle of the edge Wjbj (resp. Wj + ibj), then 

f TO + 1 (0) = (-l)^^^,^)^!.^))^^!^^^). ( 5 ) 

The cycle C corresponds to a face of the graph R*. Let c be the circumcenter of this face, and 
let Tj (resp. rj) be the angle of the rhombus R(wjbj) (resp. R{wj + ibj)) at the vertex c. Then 

Tj = ir- 28j, and rj = vr - 29'j. Since Ej=i( T i + T j) = we deduce £*L x (0j + 9 j) = ~ x )- 
Hence, 

e <Ej=i(«i+^) = _(_!)*. ( 6) 
Moreover I (wj+1 , bj) = 1 - I (6 .,^. +l) , and (-l) 1 "^^) = (-l^i+i)- 1 , so 
(-_l)E)=i(i(» i ,6 j )+i(» j+1 ,b j )) _ ^_^^Ej=i( I (» ;( . J {. i )+ I[ Ci'i^i+i)) -fc _ 



Note that Y2j=i(\wj,bj) + I(6j,«>j + i)) i s the number of co-oriented edges encountered when traveling 
cclw around the cycle C. Since the orientation of the edges of R* is clockwise odd, it is also 
counterclockwise odd, and so this number is odd. This implies 

(_l)Ej=i(i(» J ,i> J )+i(- j+1 ,b,)) _ ( 7 ) 

The proof is completed by combining equations (J5J) , © and (JJJ) . □ 

As in the complex case, a real angle in [9q — ir + A, 6$ + ir — A] can be assigned to each pole of f wb , 
for some small A > 0; where 6q is the real angle assigned to the vector b — w. 

Lemma 4.8 There exists a unique K _1 satisfying the above two properties, and K _1 is given by: 

K" 1 (b, W ) = J— I f wb ( z ) log z dz, (8) 



4vr 2 i J c 

where C is a closed contour surrounding cclw the part of the circle {e lS \ 6 G [6q — 7r + A, ^o + 7r — A]}, 
which contains all the poles off w b, and with the origin in its exterior. 

Proof: Let F(b,w) be the right hand side of ©• Fix a vertex wq G W, and let us prove that 
SfeGB K(wo, b)F(b, w) = 6 Wo (w). Denote by bi, . . . , bk the black neighbors of wq. Using Remark 14.61 
we obtain for every j, 

K(w , bj) = f Wobj (0)K(w , bj), 
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Moreover Vu> £ W, V6 £ -B, we have f«,b(0) = f«,b(0) 1 = ffc«;(0), so that f^^. (0)f TO &,. (0) = 
fw bj(0)fbjw(ty = f«;o«>(0)- Hence, using Theorem 13. 11 we obtain for every j, 

K(w ,bj)F(bj,w) = f WoW (Q)K(w ,bj)K-~ 1 (bj,w). 

Since K(wo,b) = when wq and b are not adjacent, and since K~ l is the inverse Dirac operator, 
we obtain 

k 

b£B j=l 

k 

= f™ oU ,(0)^if(wo,&j)^~ 1 (&j,?«) = f WoW (0)5 Wo (w) =S Wo (w). 

3 = 1 

Uniqueness of K _1 follows from the uniqueness of K . □ 



4.3 Convergence of the Boltzmann measure on the torus 

Let R be a rhombus tiling of the plane, and R be the corresponding rhombus-with-diagonals tiling. 
Suppose that the critical weight function is assigned to edges of R* , and denote by K the Dirac 
operator indexed by the vertices of R*. 

Consider a subset of edges e\ = wibi, . . . , = w^bk of R*, and let P be a finite simply connected 
sub-graph of R such that P* contains these edges. By Proposition 14.11 there exists a periodic 
rhombus tiling of the plane S that contains P. Let S be the corresponding rhombus-with-diagonals 
tiling, and assign the critical weight function to edges of S*. Denote by A the lattice which acts 
periodically on S, and suppose that the dual graph S* of the toroidal graph S n = S/nA is bipartite 
(this is possible by eventually replacing A by 2A). Denote by /i^ be the Boltzmann measure on 
dimer configurations A4(S^) of S*. Then we have, 



Proposition 4.9 

lim n^(ei, ... ,e k ) = [T\ K(wi,bi)) det K^fawj). 

Proof: Let us first define an orientation of the edges of S*, and the four Kasteleyn matrices 
K™, . . . , K4 of the graph 5*. Consider the graph 5J, then it is a bipartite graph on the torus. Fix a 
reference matching Mq of S^. For every other perfect matching M of <5*, consider the superposition 
M U M of M and M , then M U M consists of doubled edges and cycles. Let us define four parity 
classes for perfect matchings M of S^: (e,e) consists of perfect matchings M for which cycles of 
MDMq circle the torus an even number of times horizontally and vertically; (e,o) consists of perfect 
matchings M, for which cycles of M U Mq circle the torus an even number of times horizontally, 
and an odd number of times vertically; (o,e) and (0,0) are defined in a similar way. By Tesler |17j . 
one can construct an orientation of the edges of S*, so that the corresponding weighted adjacency 
matrix K\ has the following property: perfect matchings which belong to the same parity class have 
the same sign in the expansion of the determinant of K}. By an appropriate choice of sign, we can 
make the (e,e) class have the plus sign in det K\, and the other three have minus sign. Consider a 
horizontal and a vertical cycle of S\. Then define K3, (resp. K3) to be the matrix K\ where the sign 
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of the coefficients corresponding to edges crossing the horizontal (resp. vertical) cycle is reversed; 
and define K4 to be the matrix K\ where the sign of the coefficients corresponding to edges crossing 
both cycles are reversed. 

The orientation of the edges of §1 defines a periodic orientation of the graph S*. For every re, 
consider the graph 5* and the four matrices K™, , K3, defined as above. These matrices are 
called the Kasteleyn matrices of the graph S*. 

The orientation defined on the edges of the graph S* is a clockwise odd orientation. Let K5 be 
the real Dirac operator indexed by the vertices of S* corresponding to this clockwise orientation, 
and let Kg be the Dirac operator indexed by the vertices of S*. Then Proposition 14.91 is a direct 
consequence of Lemmas 14.101 |4~TH 14.151 below. 

Lemma 4.10 

lim //f(ei, ...,e k )= TT K 3 (wi, b t ) det K" 1 ^,^). (9) 

Proof: The toroidal partition function Z{S* L , v) is defined to be the weighted sum (weighted 
by the function v) of dimer configurations of the graph <S*. Then, by Tesler J3| (it is a generalization 
of a theorem of Kasteleyn 0), we have 

Theorem 4.11 0Q11 

Z(S*, v) = ^(-det K? + det + det K™ + det KJ). 

Kenyon gives the following theorem for the Boltzmann measure 
Theorem 4.12 7. /^(ei, . . . , eu) is equal to (jj^ =1 Kg(wi, bi) \ times 

This part of the argument can be found in 7 . Equation (jlUj) is a weighted average of the four 
quantities deti<jj<fc ((K™) _1 (6j, Wj)) , with weights 

I det K"/Z(5*, u). These weights are all in the interval (—1,1) since, for every I = 1,...,4, 
2Z(S^,u) > |detK™|. Indeed, Z(S*,u) counts the weighted sum of dimer configurations of S*, 
whereas | det K™| counts some configurations with negative sign. Moreover, by Theorem 14. 1 1 1 these 
weights sum to 1, so that the weighted average converges to the same value as each det ((K™) _1 Wj 

\<i,j<k 

Denote by B$ (resp. Ws) the set of black (resp. white) vertices of S* . Let us prove that for every 
£ = 1, ... ,4, and for every w € Ws, b € B$, (K™)~ 1 (6, w) converges to K^ 1 (6, w) on a subsequence 
of n's. The following theorem of |Hj gives the convergence on a subsequence of re's of the inverse 
Kasteleyn matrices of the graph S*. 

Theorem 4.13 (HI For every w G Ws, b G Bs, £ = 1, . . . , 4, 

r / (ien\-i(h ^ 1 f Qb,w{z,u)u x Z y dzdu 

hm n _ +00 (K7) (b,w) = — -2 / — , 11) 

(27T) 2 7 s i xS i P{z,u) z u 

where Qb, w and P are polynomials (Qb w only depends on the equivalence class of w and b), and 
x (resp. y) is the horizontal (resp. vertical) translation from the fundamental domain of b to the 
fundamental domain of w. 
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Denote by F(b,w) the right hand side of pip. In [S], it is proved that F(6, w) converges to as 
\b — w\ — » oo, as long as the dimer model is not in its frozen phase. Moreover, it is proved in ^U] 
that dimer models on isoradial graphs are never in their frozen phase when the weight function is 
the critical one. Hence, we deduce that F(6, w) converges to as \b — w\ — » oo. 
Let us prove that for every b G Bs,w & W$, F(b,w) = K^ 1 (6, w). Consider W\,W2 € W$, and 
denote by bi, . . . , the neighbors of w\. Assume n is large enough so that the graph S* contains 
W\, W2, b\, . . . , 6^, and so that the edges w\bj do not cross the horizontal and vertical cycle of S n . 
Then, for every t = 1, . . . , 4, 

k 

52 ^(w l ,b)(KrHb,w 2 ) = Y^K{^){K)-\b 3 , W 2) = S W1W2 . 

b£B s 3 = 1 

Moreover, K™(wi,bj) = Ks(wi,bj), so that taking the limit on a subsequence of n's, and using 
Theorem 14.131 we obtain 

k 

52^s(wi,bj)f{bj,w 2 ) =6 WlW2 . 
j=i 

This is true for all w\,W2 € Ws- Moreover, lim^y,^^ F(6, w) = 0, so that by definition of the 
inverse real Dirac operator, and by the existence and uniqueness Lemma 14.81 we deduce that for 
all beB s ,we W s , F(b,w) = K^ 1 (&,«;). 

Hence, /%(ei, . . . , e&) converges to the right hand side of @ on a subsequence of n's. By Sheffield's 
Theorem |16| . this is the unique limit of the Boltzmann measures u^, so that we have convergence 
for every n. □ 



Lemma 4.14 

Proof: By definition of the determinant, the left hand side of p2j) is equal to 



53 s S no " I II K s(w i ,b l ) K s 1 (fei,«) (T(1) ) . . . K s 1 (6 fc ,u) CT(fc) ), 

ctG5„ \i=l / 

where S n is the set of permutations of n elements. A permutation a G S n can be written as a 
product of disjoint cycles, so let us treat the case of each cycle separately. Refer to Section 14.2.21 
for the definition of the function f wx . 

• Suppose that in the product there is a 1-cycle, that is a point j such that a(j) = j. Then, using 
Remark 14.61 and Lemma 14.81 we obtain 



Ks(wj,bj) = f Wjb .(0)K s (wj,bj), 
K^ib^wj) = J^b~W)Ks l (b ]lWj ). 



Moreover, f^&^O) = ^^(0) 1 , hence 

^s{wj,bj)K^ l (bj,Wj) = Ksiwjib^K^^^Wj). (13) 
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• Suppose that in the product there is an ^-cycle, with £ ^ 1. To simplify notations, let us assume 
<t(1) = 2, . . . , a(£) = 1, and let us prove the following (indices are written cyclically, i.e. £ + 1 = 1), 

t I 

H K s (w j ,b j )Kg 1 (b j ,w j+1 ) = Y[K s (w J ,b J )K s 1 (b j ,w j+1 ). (14) 

3=1 3=1 

Again, using Remark 14.61 and Lemma 14.81 we obtain 

e £ 

3=1 3=1 

Using the definition of the function f wx , and the fact that it is well defined yields 

n ^bj^f^+ib^o)- 1 = n ^(0)^+1(0) = n f w+i(°) = f«. lwl (o) = 1. 
j—i j =1 j—i 

This proves equation Q14|h Combining equations ((T3*|) . Q14|). and the fact that every permutation 
is a product of cycles, we obtain Lemma 14.141 □ 



Lemma 4.15 

' T\K s (wi,bi)) det fifg 1 (6 < ,«;j)) = ( TT ] det (if-^bj, ^-)) 

2- J l<i, j<k \ / \ i / KijXfc 

Proof: Since P is simply connected, for every i,j = l,...,k, it contains a path of R from 
u>j to b{. Moreover by Theorem 13.11 the coefficient of the inverse Dirac operator correspond- 
ing to bi,Wj only depends on such a path. Hence Kg (b{,Wj) = K~ l {bi,Wj). We also have 
Vi = 1, . . . , k, Kg(wi, bi) = K(wi, so that we deduce Lemma T4. 151 □ 

□ 



4.4 Proof of Theorem CO 

The edges of the graph R* form a countable set. For every i 6 N, define fi : A4(R*) — > {0, 1} by 



MM) 



1 if the edge ej belongs to M, 
else. 



Fix k £ N, and a fc-tuple (si, . . . , s k ) of distinct elements of N. Let H € £>{0, l} fc , where £>{0, l} k 
denotes the Borel cr-field of {0, l} fc , and define a cylinder of rank k by 

A {si _ Sk) (H) = {M e M(R*)\(f Sl (M), f Sk (M)) e H}. 

Then Af ai ^ a \(H) can be written as a disjoint union of cylinder sets, 

m 

A( su ..., Sh )(H) = \J{e til ,...,et u .}, 

8=1 
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(recall that for every i, {et a , . . . , et u , } denotes the set of dimer configurations of R* containing the 
edges e Ul e tu .). Define 




Let P be a finite simply connected sub-graph of R such that, for every i = 1, ... ,m, P* contains 
the edges et il , . . . , et u , ■ Let S be a periodic rhombus tiling of the plane that contains P (given by 
Proposition 14. lj) . Then, by Proposition 14.91 [ir Su 8k )(H) = lim n _,. 00 /i„(A( Sl) ___ !Sfe )(iJ)). From this 
we deduce that for every k, and for every fc-tuple (sx, . . . , sjt), M( Sll ... )S& ) is a probability measure on 
B{0, l} k . Moreover, we deduce that the system of measures {/i( si) .„ )8fc ) : (si, • • • , Sfc) is a /c-tuple of 
distinct elements of N} satisfy Kolmogorov's two consistency conditions. Applying Kolmogorov's 
extension theorem, we obtain the existence of a unique measure fj, R , which satisfies (|2*j). 

Using the fact that the measure ji R of a cylinder set is the limit of Boltzmann measures, we deduce 
that the measure n R is a Gibbs measure in the sense given in the introduction. 

Assume moreover that the graph R* is doubly periodic. Then, for every cylinder set {ex, . . . , e&} 
of R*, /J, R (ei, . . . , efe) = limfi-Kx, fj, R (ex, . . . , e&). Moreover by ^Hl (see also the Boltzmann 
measures fi R converge to the minimal free energy per fundamental domain Gibbs measure which is 
unique, so that this proves Theorem 13.21 □ 



5 Gibbs measure on the set of all triangular quadri-tilings 

Recall that Q is the set of all triangular quadri-tilings up to isometry, i.e. the set of all quadri- 
tilings whose underlying tiling is a lozenge tiling of the equilateral triangular lattice T. Denote by 
A4 the set of dimer configurations corresponding to quadri-tilings of Q. In this section, we first 
define the notion of Gibbs measure on ftA; then we define a <r-algebra cr(B) on A4, and give an 
explicit expression for a Gibbs measure /i on (M,a(E)). We conjecture fi to be of minimal total 
free energy per fundamental domain among a four parameter family of ergodic Gibbs measures. 

The notion of Gibbs measure on Q is a natural extension of the notion of Gibbs measure on dimer 
configurations of a fixed graph. Assume a weight function v is assigned to quadri-tiles of triangular 
quadri-tilings of Q, then a Gibbs measure on Q is a probability measure with the following prop- 
erty. If a triangular quadri-tiling is fixed in an annular region, then quadri-tilings inside and outside 
of the annulus are independent. Moreover, the probability of any interior triangular quadri-tiling is 
proportional to the product of the weights of the quadri-tiles. Using the bijection between Q and 
M, this yields the definition of a Gibbs measure on Ai. 

Define C to be the set of lozenge-with-diagonals tilings of the plane, up to isometry. Define C* to 
be the graph (which is not planar) obtained by superposing the dual graphs L* of lozenge-with- 
diagonals tilings L € C. Although some edges of C* have length 0, we think of them as edges of 
the one skeleton of the graph, so that to every edge of C* there corresponds a unique quadri-tile 
in a lozenge-with-diagonals tiling of C. Let e be an edge of £*, and let q e be the corresponding 
quadri-tile, then q e is made of two adjacent right triangles. If the two triangles share the hypotenuse 
edge, they belong to two adjacent lozenges; else if they share a leg, they belong to the same lozenge. 
Let us call these lozenge(s) the lozenge(s) associated to the edge e, and denote it/them by l e 
(that is l e consists of either one or two lozenges). Let k e be the edge(s) of T* corresponding to the 
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lozenge(s) l e . 

Let e±, . . . , ef. be a subset of edges of £*, and define the cylinder set {ei, . . . , e&} of C* to be the 
set of dimer configurations of hA which contain these edges. Let us call connected cylinder any 
cylinder of C* which has the property that the lozenge(s) associated to its edges form a connected 
path. Then every cylinder of C* can be expressed as a disjoint union of connected cylinders. Con- 
sider B the field consisting of the empty set and of the finite disjoint unions of connected cylinders. 
Denote by cr(B) the a-field generated by B. 

Let fi be the Gibbs measure on dimer configurations of the honeycomb lattice A4(T*) given in 
[7]. As we have noted before, Theorem 13.21 is true for all doubly periodic isoradial graphs with 
critical weights on their edges. Then, the measure ^ coincides with the minimal free energy per 
fundamental domain Gibbs measure given by Theorem 13.21 when the doubly periodic isoradial 
graph is the honeycomb lattice. As a corollary to Theorem 13.21 we obtain: 

Corollary 5.1 There is a probability measure \x on (Ai,a(B)) such that for every connected cylin- 
der {ex, . .. ,e k } of £*, 

H(ex,...,e k ) =/x L (ei,...,e fc )/i T (k ei ,...,k e J, (15) 

where L is the lozenge- with- diagonals tiling corresponding to any lozenge tiling L which contains 
the lozenges l ei , . . . , l efc . Moreover fj, is a Gibbs measure on M., where the critical weight function 
is assigned to quadri-tiles. 

Proof: Expression (|15fl is well defined, i.e. independent of the lozenge tiling L which contains the 
lozenges l ei , . . . , l efc . Indeed, by definition of a connected cylinder set, the lozenges associated to the 
edges ei, . . . , e k form a connected path of lozenges, say 7. Let L be a lozenge tiling that contains 
7, and denote by Kl the complex Dirac operator indexed by the vertices of the graph L*. Then 
K^ibi^Wj) is independent of the lozenge tiling L which contains 7, indeed K^fa, Wj) only depends 
on an edge-path of R from Wj to bi, and since L contains 7 which is connected, we can choose the 
edge-path to be the same for all such lozenge-with-diagonals tilings L. We then use the fact that 
/i T and fi L are probability measures to prove the two conditions of Kolmogorov's extension theorem. 

The measure is a Gibbs measure on M. as a consequence of the fact that ji L and fjJ are Gibbs 
measure on A4(L*) and A4(T*) respectively. □ 

Assume that a weight function is assigned to quadri-tiles of quadri-tilings of Q. Denote by A the 
lattice which acts periodically on T, and by T n = T/raA, moreover suppose that T* is bipartite (this 
is possible by eventually replacing A by 2A). Define Q n to be the set of triangular quadri-tilings 
whose underlying tiling is a lozenge tiling of T ra . Denote by fi n the Boltzmann measure on Q n ; 
that is the probability of having a given subset of quadri-tiles in a quadri-tilings of Q n chosen with 
respect to fi n is proportional to the product of the weights of the quadri-tiles. We make the first 
following conjecture. 

Conjecture 1 Suppose the critical weight function is assigned to quadri-tiles, then the Gibbs mea- 
sure of Gorollarv \5.1\ is the limit of the Boltzmann measures \x n . 

Now fix (s,t,p,q) S M 4 , and let Q^ s '*' p ''^ be the subset of quadri-tilings of Q n whose first height 
change is (l/isj, L n ^J) and second height change is ([np\, [nq\). Assuming that Q^ ,t,p,q ^ is non 
empty, let (j£ ,t,p,q ^ be the conditional measure induced by fi n on d sAM) . Denote by the 



25 



weighted sum of quadri-tilings of Qn ,t,p,q \ and define the total free energy per fundamental 
domain a(s,t,p,q) by: 

a(s,t,p,q) = - lim -I log 

n— >oo n 

Then the following conjecture is inspired by a result of ^B], see also PJ. 

Conjecture 2 For eac/i (s,t,p,q) for which Q^ ,t,p,q ' is non empty for n sufficiently large, fin' t ' p ' qS> 
converges as n — > oo to an ergodic Gibbs measure fj,( s > t ' p ' q > of slope (s,t,p,q). Furthermore fj> n itself 
converges to ^( s o,*o,po,<?o) w h ere ( S0; i ,p , q ) is the limit of the slopes of fi n . If (sq, to,po, qo) lies in 
the interior of the set (s,t,p,q) for which Q£> t > p > q > is non-empty for n sufficiently large, then every 
ergodic Gibbs measure of slope (s,t,p,q) is of the form ^ S ^'P^) for some (s,t,p,q) as above; that is 
^(f,t,p,q) j fi un iq Ue ergodic Gibbs measure of slope (s,t,p,q). Moreover, the measure ^( s o>*o.po,9o) 
is the unique one which has minimal total free energy per fundamental domain. 

As a consequence, when the critical weight function is assigned to quadri-tiles, we conjecture the 
minimal total free energy per fundamental domain Gibbs measure to be given by the explicit 
expression of Corollary 15.11 

6 Asymptotics in the case of quadri-tilings 

Section fn. II aims at giving a precise statement of Theorem ll.2l of the introduction (see Theorem 16, II 
below) , and Section 16.21 gives consequences of Theorem 16.11 for the Gibbs measure fj, R of Theorem 
K-}.2l and /i of Corollary 15.11 

6.1 Asymptotics of the inverse Dirac operator 

Refer to Figure ITTI for the following notations. Let £ 2 be two disjoint side-length two rhombi 
in the plane, and let £\, £2 be the corresponding rhombi with-diagonals. Assume l\ and £2 have a 
fixed black and white bipartite coloring of their faces. Let r% and 7*2 be the dual graphs of i\ and 
£2 ( r i and T2 are rectangles), with the corresponding bipartite coloring of the vertices. Let robea 
white vertex of 7*1, and b a black vertex of 7*2, then w (resp. b) belongs to a boundary edge e\ of i\ 
(resp. e2 of £2)- By Lemma l2~2l to the bipartite coloring of the faces of £\ and £2, there corresponds 
a bipartite coloring of the vertices of £[ and £' 2 . Let x\ (resp. X2) be the black vertex of the edge 
ei (resp. e%). Orient the edge wx\ from w to x\, and let e tSl be the corresponding vector. Orient 
the edge X2b from X2 to 6, and let e 1 ® 2 be the corresponding vector. Assume l\ and £' 2 belong to a 
rhombus tiling of the plane R. Moreover, suppose that the bipartite coloring of the vertices of R* 
is compatible with the bipartite coloring of the vertices of r\ and 7"2 . 





Figure 11: Rhombi with diagonals £±, £2 and their dual graphs r%, r%- 
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Then we have the following asymptotics for the inverse Dirac operator K indexed by the vertices 
of R*. Refer to the introduction for comments about Theorem 16.11 

Theorem 6.1 As \b — w\ — ► cxd, K~ 1 (b,w) is equal to 

J_ ( 1 , e -'< fl l+ e 2) \ J_ / e 2, 9l+e 2, 92 e -i(39 1+ fl 2 ) +e -i(« 1+ 38 2 ) \ / 1 \ 

2tt ^b— tu "T b-w J ~*~ 2tt \ (f>-io) a h (b—ffi) 3 J V I 6 -™! 3 / ' 

where Q\ and 62 are defined above. 

Proof: Let us define an edge-path from w to b in R (the set of rhombi associated to the edges 
of R*). Consider the bipartite coloring of the vertices of R (given by Lemma l2.2j) associated to the 
bipartite coloring of the vertices of R* . We define the graph N as follows. Vertices of N are black 
vertices of R, and two vertices of N are connected by an edge if they belong to the same rhombus 
in R. The graph N is connected because R is. Each face of N is inscribable in a circle of radius two. 
The circumcenter of a face of N is the intersection of the rhombi in R, to which the edges on the 
boundary cycle of the face belong. Thus the circumcenter is in the closure of the face, and so faces 
of N are convex. Note that the vertices x\ and X2 are vertices of the graph N. 

Denote by (x, y) the line segment from a vertex x to a vertex y of N. An edge uv of N is called 
a forward-edge for the segment (x, y) if < v — u, y — x > > 0. An edge-path v±, . . . ,Vk of N is 
called a forward-path for the segment (x,y), if all the edges iro+i are forward-edges for (x,y). 
Similarly to what has been done in [H], let us define a forward-path of N for the segment (x\,X2), 
from x\ to X2 (see Figure EJ- Let Fi, ...,F( be the faces of N whose interior intersect (xi,x>z) 
(if some edge of N lies exactly on (xi,X2), perturb the segment (xi,X2) slightly, using instead a 
segment (x\ -\-E\,X2 + £2) for two generic infinitesimal translations (ex, £2))- Note that the number 
of such faces is finite because the rhombus tiling of the plane R has only finitely many different 
rhombi. Then for j = 1, . . . ,£ — 1, Fj n Fj+i is an edge e^+i of N crossing {x\,X2)- Set v\ = x\, 
V£ = X2, and for j = 1, ...,£ — 2, define Vj + \ to be the vertex of ej + \ such that the edge e^+i 
oriented towards Vj + \ is a forward-edge for {x\,x<i). Then, for j = 1, . . . , t — 1, the vertices Vj and 
Wj-4.1 belong to the face Fj. Take an edge-path from Vj to Vj+i on the boundary cycle of Fj, such 
that it is a forward-path for (xi,xz). Such a path exists because faces of N are convex. Thus, we 
have built a forward-path of N for (x\, X2), from x\ to X2- Denote by u\ = x\, u^, ■ ■ ■ , u^i, = X2 
the vertices of this path. 




Figure 12: forward-path from x\ to X2 for the segment {x\,X2j. 
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Let us now define an edge-path of R from w to b. Note that the edges wx\ and x 2 b are edges of R. 
For j = 1, . . . , k — 1, define the following edge-path of i? from Uj to %+i (see Figure EU- Remember 



that UjUj + \ is the diagonal of a rhombus of R, say 



Let fj 



be the dual graph of 



j- 



Let Uj 



be 



the black vertex in fj adjacent to uj, let u* be the crossing of the diagonals of and let iq be the 



white vertex in fj adjacent to Uj+i. Then the path iij,Uj,u 2 



_i is an edge-path of R. Thus 



w,xi = Ul,U{,Uf,Ul,U2, ■ ■ ■ ,u k -i,ul 



v u k- 



l> 



u k 



X2 , b is an edge-path of R, from u) to 6. 



Orient the edges in the path towards the black vertices of R* , and away from the white vertices of 
R*. 



u 


1 / 
j 




3 






*-s 





















Figure 13: Edge-path of R from to Uj + \. 

Let e 1 ' 3 ' , e^i , e a i , e lQ ? be the vectors corresponding respectively to the edges Ujuj, i$Uj+i, ujuj, 
u^u 2 . Without loss of generality, suppose that x 2 —x\ is real and positive. Then for j = 1, . . . , k — 1, 
and £ = 1,2, we have: 



at I < ~ Uj,X 2 ~ Xi > 

cos m - cos oq = j — j . 

2|X2 — X\\ 

Since u±, . . . ,u k is a forward-path for (xi,x 2 ), this quantity is positive, thus cos/?j > cosaj. 
Moreover, since there is only a finite number of different rhombi in R, k = 0(\b — w\). For the same 
reason, there is a finite number of angles ffi, and they are all in [— tt + A,ir — A], for some small 
A > (in the general case where the angle of the vector x 2 — X\ is Oq, the angles m would be in 

the interval [9o — tt + A, 6>o + tt — A]). Thus by Theorem 4.3 of 8. , we have that K _1 (6, w) is equal 
to: 

1 ^ ■ 1 ' 6 ■ * ^-O('). (16) 



27r \b — w ' b — wj 2tt \(b — w) 3 (b — w) 3 J \\b — w\ 3 
k—1 2 k—1 2 

where 7 = e -*(*i+ft») J] f[ e i ^ +a '\ and £ 2 = e 2 ^ + e 2 ^ + Y,Y1 
j'=l £=1 j=l i=l 

Note that for j = 1, . . . , A; — 1, we have a 2 = (/?] + tt) mod[27r], and (3 2 = (a) + tt) mod[27r], thus: 



2 



n 

2 

53 e 2i ^ - e 2ia i = e 2i(3 i - e 2ia ) + e 2l K+-) - e 2l( ^ +n) = 0. 
i=i 

Therefore 7 = e - i( - 9l+d2 \ & = e 2idl + e™ 2 , which proves the theorem. □ 
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6.2 Asymptotics of the Gibbs measures on quadri-tilings 

Let R be a rhombus tiling of the plane, and R be the corresponding rhombus-with-diagonals tiling. 
Consider a subset of edges e% = itfi&i, . . . ,e k = w k b k of R* , and recall that fi R is the Gibbs measure 
on Ai(R*) given by Theorem 13.21 

Corollary 6.2 When Vj 7^ i, \wj — b%\ — > 00, then up to the second order term, fx R (ei, . . . , e k ) only 
depends on the rhombi of R to which the vertices bi and Wj belong, and else is independent of the 
structure of the graph R. 

Proof: This is a consequence of the explicit formula for fi R (e\, . . . ,e k ) of Theorem 13.21 and of 
the asymptotic formula for the inverse Dirac operator of Theorem 16.11 □ 

Recall that C* is the non-planar graph obtained by superposing duals of lozenge-with-diagonals 
tilings of C. Let e\ = W\b\, . .. ,e k = w k b k be a subset of edges of C*. Define C E to be the set 
of lozenge-with-diagonals tilings of the plane that contain the lozenges associated to the edges 
ei, ... ,e k . 

Corollary 6.3 When Vj 7^ i, \wj — 6j| — > 00, then up the second order term, fA L (ei, . . . , e k ) is 
independent of L S C E . 

Proof: As in Section we choose an embedding of C* so that every edge of C* uniquely deter- 
mines the lozenge(s) it belongs to. Corollary 16.31 is then a restatement of Corollary 16.21 □ 

Recall that \i is the Gibbs on triangular quadri-tilings given by Corollary 15. II Let \ ei , . . . ,\ eh be the 
lozenges associated to the edges ei, . . . , of £*, and let k ei , . . . , k efc be the corresponding edges of 
T*. 

Corollary 6.4 When Vj 7^ i, — 6j| —> 00, and for every L € C E , we have 



Proof: This is a consequence of the explicit formula for /i(ei, . . . , e^) of Corollary 15.11 and of 
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